Quantum-state steering in optomechanical devices 
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We show that optomechanical systems in the quantum regime can be used to demonstrate EPR-type quantum 
entanglement between the optical field and the mechanical oscillator, via quantum-state steering. Namely, the 
conditional quantum state of the mechanical oscillator can be steered into different quantum states depending 
the choice made on which quadrature of the out-going field is to be measured via homodyne detection. More 
specifically, if quantum radiation pressure force dominates over thermal force, the oscillator's quantum state is 
steerable with a photodetection efficiency as low as 50%, approaching the ideal limit shown by Wiseman and 
Gambetta [Phys. Rev. Lett. 108, 220402 (2012)]. We also show that requirement for steerability is the same as 
those for achieving sub-Heisenberg state tomography using the same experimental setup. 



Introduction. — Optomechanical devices are now approach- 
ing the quantum regime and can therefore provide a platform 
for investigating quantum behaviors of macroscopic mechan- 
ical degrees of freedom [Ij. A topic extensively discussed 
is how to demonstrate quantum entanglement in such sys- 
tems GrIH). In this paper, we focus on an interesting aspect of 
bipartite quantum entanglement — the ability of modifying the 
quantum state of one party by making different measurements 
on the other. This is essence of the Gedankenexperiment of 
Einstein, Podolsky, and Rosen (EPR) |9|, and has been rigor- 
ously formulated by Wiseman et al. ifToHL?! as quantum steer- 
ability. More recently, Wiseman has shown that steerability 
can be demonstrated by showing detector-dependent stochas- 
tic evolution of a two-level atom coupled to an optical field 
which in turn is measured continuously 1T31 . 

In the context of linear quantum systems with Gaussian 
states (appropriate for optomechanics experiments), steering 
can be understood as follows. According to quantum mechan- 
ics, position and momentum of a mechanical oscillator satisfy 
the Heisenberg uncertainty principle, which reads: 



AX^ AX fc > |sin(0i-02)|, V0i,0 2 



(1) 



where = (x/Ax q ) sin0 + (/3/A/^) cos 0, with Ax q and 
Ap q zero-point uncertainties in position and momentum, are 
quadratures of the mechanical oscillator. Let us assume a 
mechanical oscillator is interacting and establishing entangle- 
ment with a continuous optical field, and then the field gets 
continuously measured via time-dependent homodyne detec- 
tion with 0(f) being the optical quadrature measured at time 
t. Suppose the measurement lasts from — T up to 0, the final 

conditional state of the oscillator, written as (0)), will de- 
pend on how we make the homodyne detection due to entan- 
glement. For two different measurement strategies with 6\ (t) 
and 02(f), respectively, in general we have two different final 

conditional states: |y/m'(0)) 7^ | y/,',f 2 (0)). If the quadratures 
are properly chosen, we may have, as illustrated in Fig.[T] 
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where AX^ = ((% - (X^.)) 2 ) 1 ^ with (•) = (y/,7| • In 
other words, if in the first strategy, the observer tries to predict 



quadrature X^ of the mechanical oscillator, while in the sec- 
ond strategy, the observer tries to predict X<^, then the two pre- 
dictions have an error product that is lower than Heisenberg 
Uncertainty. This is the essential feature of the EPR paradox; 
the ability of an experimental setup in creating such a feature 
is called quantum steerability II 1 044 1 21 . 

In ideal linear quantum measurement processes, both con- 
ditional states will be pure, and for any pairs of distinctive 
0i and 02, inequality (|2) will almost always exist for some 
set of 0i and 02 — although this idealized steerability may be 
influenced by practical imperfections, such as thermal noise. 
In view of Eq.|2|, we introduce a figure of merit to quantify 
steerability, 



5? = — min 

01 ,02,9i ,62 



|sin(0i -0 2 )| 



(3) 



with minimum obtained by comparing all possible sets 
of {0i, 0i (?), 02, 2 (f)} {t G 0]) — an optimal time- 
dependent homodyne detection is needed to achieve the lower 
bound. The quantum state is steerable when 5? > 0, which 
will be proved to be equivalent to the formal criterion obtained 
by Wiseman for Gaussian entangled states ifTUl . As we will 
show in the discussion that follows, for linear optomechani- 
cal devices, when the quantum radiation pressure dominates 
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FIG. 1: (color online) Two different quantum states (projection of 
their Wigner functions on phase space) of a mechanical oscilla- 
tor, conditional on two different strategies for measuring the optical 
quadrature: one at 8\ (t) (left) and the other at 82(f) (right). 
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strongly over thermal fluctuations, steerability only depends 
on the photodetector efficiency 7} of time-dependent homo- 
dyne detections: 



2 [lnT? 



-ln(l -!?)], 



(4) 



which will be positive as long as r\ > 50%, which coincides 
with the ideal limit shown by Wiseman and Gambetta lfl3l . 
Interestingly, such quantum steerability is intimately related 
to the state tomography accuracy in the protocol suggested by 
Miao et al. [14], in which an optimal time-dependent homo- 
dyne detection scheme is used to probe the quantum state of a 
mechanical oscillator with Gaussian-distributed joint position 
and momentum error less than Heisenberg uncertainty. More 
explicitly, we will show, for the same optomechanical device, 



In IVdetVy/Tz 



(5) 



where V,, is the covariance matrix for the tomograph error. 

Optomechanical dynamics. — We start by considering dy- 
namics of linear optomechanical device which has been ex- 
tensively studied ll8l [131117 1 : its linearized Hamiltonian reads: 

3$ =p 2 /(2m)+malx 2 /2 + j£ Km +hAa Jf a + hgx(a Jf +a) 
+ ihy/K(a ext a^ — fle Xt fl) • (6) 

Here co m is the mechanical resonant frequency; A = co c . — ©/ 
is cavity detuning, i.e., difference between the cavity reso- 
nant frequency co c and the laser frequency cor, ^k„, summa- 
rizes the fluctuation-dissipation mechanism for the mechan- 
ical oscillator; the fifth term is the optomechanical interac- 
tion term with g = aco c /L quantifying the coupling strength, 
a the steady-state amplitude of the cavity mode and L the cav- 
ity length; the last term describes the coupling between the 
cavity mode and external continuous optical field a ext with 
[d ext (t), a ext (f')] = S(t — t') in the Makovian limit and K is the 
coupling rate which is also the cavity bandwidth. 

From Hamiltonian |6]), one can obtain the following set of 
linear Heisenberg equations of motion 

m[i(f ) + K m k{t) + 0), 2 „i(f)] = Prp(t) + A(0 . (7) 
&(t) + (k/2 + iA)d(t) = -igx(t) + VKd in (t) , (8) 



and the input-output relation 



flout (0 = -&m(t) + Vxd(t) , 



(9) 



where d in = a ext (f-) (in-going) and a out = a ext (f+) (out-going) 
are input and output operators in the standard input-output for- 
malism 1 18 1, and F^, = —hg(a + eft) is the quantum radiation 
pressure force and is the thermal fluctuation force associ- 
ated with the mechanical damping. Output field a out can be 
measured with a homodyne detection scheme, from which we 
can infer the mechanical motion and thus the quantum state of 
the oscillator. By adjusting the local oscillator phase, one can 
measure any ^-quadrature: bg = b\ sin 9 + bi cos 9, which is 
a linear combination of the output amplitude quadrature b\ = 
(flout + flout) and P hase quadrature b 2 = 4^. (flout - flJut)- 



After incorporating non-unity photodetection efficiency r\ , the 
measurement output at time t is given by 



9e (0 = y/rf [*i (0 sin 9 + bi(t) cos 9} + y/l-rj ne (?) , (10) 

where fig is the vacuum noise associated with the photode- 
tection loss and is uncorrelated with «;„. Note that 9 can 
be a function of time, when the local oscillator phase is ad- 
justed during the measurement; in this way a different optical 
quadrature (but only one) is measured at each moment of time. 

Conditional quantum state. — Suppose we perform out ho- 
modyne detection during — T < t < 0, obtaining a data string: 

ye = {ye(-r),yg(-T + dt), ••• ,y e {-dt),y e {0)} (11) 

with dt = t/(N — 1) being the time increment and the num- 
ber of data points. We can then infer the quantum state of 
mechanical oscillator at t = conditional on these measure- 
ment data, obtaining the so-called conditional quantum state. 
The standard way to obtain the conditional state is to use data 
to drive the stochastic master equation 1 19-24]. Here we use a 
different approach by using the Wigner quasi-probability dis- 
tributions, deriving the conditional quantum state in a way 
similar to classical Bayesian statistics, in the same spirit as the 
approach applied in Ref. Il25ll : this allows us to more straight- 
forwardly treat non-Markovianity, e.g., due to a small cavity 
bandwidth and colored classical noises. Our approach takes 
the advantage of the following facts: 

[ye(t),9e(t')} - lx(0),9e(t)} = \p(fi),Se{t)]=0 (12) 

Vf,?' £ [— T,0], which is a consequence of the general fea- 
tures of linear continuous quantum measurements ll26l . We 
can therefore treat yg (t) almost as classical quantities and ig- 
nore their time-ordering in deriving the following joint Wigner 
function of the oscillator and the continuous optical field: 

W(x,y e ) = np(-t)8^(x-x)8^(yg- y g)} . (13) 

Here we have only included the marginal distribution for the 
optical quadrature y g of interests, instead of the entire opti- 
cal phase space; x = (x(0), p(Q)) and x is a c-number vector, 
similar for;y e ; p(— t) is the initial joint density matrix 



|0)(0| 



(14) 



with p* the thermal state of the oscillator and |0) the vacuum 
state for the optical field — the coherent amplitude of the laser 
has been absorbed into the optomechanical coupling constant 
g. Since jc(0) and p(Q) do not commute we have to explic- 
itly define 8^ 2 \x— x) = J d 2 %e~'*'ft~ x \ Similar to classical 
Bayesian statistics, the Wigner function for the conditional 
quantum state of the mechanical oscillator at t = reads: 



W m (x\y e ) = W(x,y e )/W(y e ) 



(15) 



Since we consider only Gaussian quantum states, the joint 
Wigner function can thus be formally written as: 



W(x,y ) =c exp 



-U*>ye)y~\*,ye) 



(16) 
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where cq is the normalization factor and superscript T denotes 
transpose. Elements of the covariance matrix V are given by 



'jk 



(XjX k ) sym = Tr[p(-T)(XjX k +X k Xj)]/2 (17) 



with X = (x,y g ). We separate components of the oscillator 
and the optical field, and rewrite the covariance matrix V as: 



V = 



'ACJ" 




_ Cq Bg 





A C T u T g 

Uf,C UllKUa 



(18) 



Here A is a 2 x 2 covariance matrix for the mechanical oscilla- 
tor position i(0) and momentum p(0); B is a 2N x 2N covari- 
ance matrix for two quadratures >*i = 9e=n/2 an( l 5*2 = 5*9=0 
of the optical field; u g — (sin 6, cos0) is a N x 2N matrix 
and sinO = diag[sin 0(— t), • • • ,sin0(O)] — a diagonal matrix 
with elements being quadrature angle at different times; C is a 
2N x 2 matrix describing the correlation between (5*1,5*2) and 
(i(0), p(Q)). Combining Eqs. ( fl"5] l and ([l8j, we obtain: 



W m (x\y e ) = — exp 



- l -(x-x\ e )\\r\x-xn 



9\T 



, (19) 



where the conditional mean x> and covariance matrix Vj„ are 



x^=C T e B g y e , V, 



— A — CnBo'Cfi 



(20) 



Note that the two rows of the 2N x 2 matrix, CgBg , which 
we shall refer to as K x (the first row) and K p (the second 
row), are also the optimal filters that predict i(0) and p(0) 
with minimum errors, ([jc(0) — K x y1] 2 ) and (\p(0) — K p y T e ] 2 ), 
respectively. The above results for the conditional mean and 
variance are formally identical to those obtained by classical 
optimal filtering. 

Quantum-state steering. — We can now understand the 
quantum-state steering from a more quantitative way. From 
Eq. ( pO) , we learned that the conditional variance of the os- 

cillator state V)„ directly depends on the optical quadrature 9 
that we choose to measure. To calculate the steerability figure 
of merit 5? [cf. Eq. Q], we need to find the time-dependent 
quadrature phase 0(t) that minimize the conditional variance 

AX^ of a given mechanical quadrature Xq = v,px with vector 
= (sin(j>/Axq, cos (j) /Ap q ). Using the fact that 

min(Axl e ) 2 =min(v^ T -^> 2 

= min(v^ T -^T-^ 2 yI) 2 (2D 

with K\ = Ks\n and K2 = Kcos 6, we obtain the minimum 

(Axl 9 ) 2 nin -v^(A-C T B- 1 C)vJ, (22) 

and 0(t k ) at t = — T + kdt is given by: 

8(t k ) = arctan [(v^C J B- l ) k /(v^C T B- l ) N+k \ . (23) 

Since (Ax], 9 )^ is in a quadratic form of va, we obtain: 

Y = -\n \2y/det\ s /h] , V s = A — C T B _1 C . (24) 



This means quantum state of the oscillator is not steerable — 
< 0, if Vi is Heisenberg limited — VdetV^ > h/2. 
Such a definition of steerability is in accord with the crite- 
rion by Wiseman et al. (T0|, more specifically, shown in their 
Eq. (17), which says that quantum state of the oscillator can- 
not be steered by the optical field, if we have 



A C T 
C B 



•i£ w ©0 o >0, 



(25) 



where £„, is the 2x2 symplectic matrix for the oscillator, and 
0„ is a null 2N x 2N matrix for the optical field. Since the 
covariance matrix B for the optical field is positive definite, 
namely, B > 0, the above condition requires that the Shur's 
complement of A be positive definite: 



A — C B C + iZ m =Y S + iZm > , 



(26) 



which is equivalent to requiring that Vj is Heisenberg limited, 
i.e., 



In [2v / det\y/zl <0 



(27) 



Continuous-time limit. — To properly describe the actual 
continuous measurement process, we take the continuous- 
time limit with dt — > 0, and we have — » °°. The matrices 
indexed by time become functions of time, while matrix prod- 
ucts involving summing over time become integrals. In partic- 
ular, the central problem of calculating K — C T B~' becomes 
solving an integral equation for K: 



J° dt'B{t,t')K(t')=C T {t). 



(28) 



More specifically, B(f , f') now becomes a 2 x 2 matrix with 
elements being the two-time correlation functions between 
optical quadratures y\ (t) and y2(t')\ C(f )'s elements are cor- 
relation functions between (y\ (t),y2(t r )) and (i(0),/3(0)). 
These correlation functions can in turn be obtained by solv- 
ing Heisenberg equations of motion [cf. Eqs. fTpO])] and 
expressing x(0), p(0), fi(t), and $2(t) in terms of a\ n (t), 
n(t) and F t h(t), for which we have (ai n (t)dj n (t')) sym = 8(t — 

0/2, (n(t)n(t')) sym - d(t-t')/2 and (F t h(f)^th(f'))sym = 
2mK m k.BT 8(t — t') given the initial state p(— t) shown in 
Eq. (14) . The above integral equation is generally difficult to 
solve analytically if z is finite. Since usually we are not in- 
terested in the transient dynamics, we can extend — T to — 00, 
which physically corresponds to waiting long enough till the 
mechanical oscillator approaches a steady state, and then start 
state preparation. In this case, Eq. ( f28| can be solved ana- 
lytical using the Wiener-Hopf method of which the detail is 
shown in the Appendix|A] 

Large cavity bandwidth and strong measurement limit. — 
One interesting scenario that allows a nice compact-form so- 
lution is when the cavity bandwidth is large and the optome- 
chanical coupling rate is strong — a strong measurement, com- 
pared with the mechanical resonant frequency co m . In this 
case, the cavity mode can be adiabatically eliminated, and the 
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mechanical resonant frequency co m ignored (for general sce- 
narios, the formalism here still applies but the analytical re- 
sults become quite complicated). Correspondingly, equations 
of motion for the oscillator simply reads [cf. Eqs. |7) and 



Steerability 



m'i(t) = F^(t) + Fth(f) = -aai(f) + F th (t) , 



(29) 



with d\ = 4| (a; n + ajj) the amplitude quadrature of the input 
field. The output amplitude quadrature y\ and phase quadra- 
ture $2 are given by [cf. Eqs. Q and <jT0]»] : 



9l (0 =y/n&\ if) + y/l -??«i(0 , 

hit) = Vn [02(0 + (a/h)x(t)] + vT" 



(30) 

^2(0, (31) 



where 02 = (fl; n — a^) is the phase quadrature of the input 
field, and additionally we have introduced an effective cou- 
pling constant a = ^8/Khg. From the above equations, we 
can easily obtain those correlation functions in the integral 
equation shown in Eq. (f28]>, solving which leads to: 



V2ri 



2?V« 2 /(CfM 1 

1 ll^J^phmja 2 



where the characteristic constant C,p is defined as: 



rj / AmK m k B T 



1/2 



Correspondingly, we obtain the steerability [cf. Eq. |24 
y = -ln(V2C F /ri). 



(32) 



(33) 



(34) 



For a strong measurement, the quantum radiation pressure 
dominates over the thermal fluctuation force, and we have 
5j? = a 2 > = 4mK m k B T, with and being the single- 
sided spectra density — twice the Fourier transform of two- 
time correlation function. This leads to £/? ps yT}(T— J ?)/2 
and ^ s» jln[T}/(l — 77)], as shown in Eq. Q. 

Connection between steering and quantum tomography. — 
Interestingly, such quantum-state steering is closely related 
to the quantum tomography protocol discussed in Ref. lfT4l . 
where an optimal time-dependent homodyne detection is pro- 
posed to minimize the error in obtaining marginal distribu- 
tions of different mechanical quadratures, from which we re- 
construct the Wigner function of the quantum state in phase 
space. More specifically, for the same optomechanical device 
discussed above, the tomography error — quantifying the dif- 
ference between the reconstructed Wigner function and the 
actual one — is given by the following covariance matrix: 



V v = 



Hf 
\/2?7 



2*y/a 2 /{£ F Hm) -1 

-1 l^^^phmja 2 



(35) 



Notice that it is almost identical to the conditional covariance 
matrix V. v shown in Eq. ( |32] i, apart from that the off-diagonal 
terms have the opposite sign. The state steering can there- 
fore be viewed as the time-reversal counterpart for state to- 
mography, as the off -diagonal term flips sign when the oscil- 
lator momentum p — > — p under t —> —t, and the condition 




Frequency 



10 20 30 40 50 

O x /fljr or 277/(1 - rj) 



FIG. 2: (color online) The displacement noise spectrum of an op- 
tomechanical device and the characteristic frequencies of the force 
noise Sip, the quantum noise Q. q , and the sensing noise Sl x referring 
to the Standard Quantum Limit (SQL) (left) and the steerability y v 
as a function of ratio £l x /Qp while keeping Q. q = i/Sl x £lp (right). 
Only the ratio between these frequencies determines the steerability, 
and the absolute values can vary depending on the chosen device. 



for achieving a sub-Heisenberg error for state tomography — 
V / deTVv < h/2, is also identical to that for steerability. 

Such a connection can be understood from the fact that for 
steering, one tries to prepare states with minimal uncertainty 
in certain quadratures (0) with data from (— «>, 0] — a filter- 
ing process, while for tomography, one tries to minimize the 
error in estimating quadratures X^, (0) with data from (0, °°) — 
a retrodiction process. Due to linearity, both the minimal un- 
certainty and the tomography error for a given quadrature X^ 
all takes the quadratic form — v V Sl v vl and 



V,=A-C T B 'C 



t ->• -f 



> V v . 



(36) 



These two covariance matrixes V ljV describe the remaining 
uncertainty in oscillator position x(0) and momentum /3(0) 
conditional on both the amplitude y\ and phase quadrature $2 
of optical field, for t < and t > 0, respectively. Note that the 
above relation shown in Eq. ([36| is exact only when the noise 
during the state preparation and the one during tomography 
are uncorrelated, as the correlation between them will break 
down the time-reversal symmetry, which happens if the cavity 
bandwidth is small and has noneligible memory time. 

Verifiable steering and the experimental requirements. — 
Not only are the steering and tomography intimately related 
to each other, but also the tomography is necessary in order to 
verify the steering in the experiment. For Gaussian states, the 
tomography error simply adds on top of the covariance matrix 
for every conditional state. We therefore define the following 
figure of merit for verifiable quantum-state steering: 



In 2 v /det[V J +V v ]//z = -ln(2fr /17). (37) 



We therefore require £p < rj/2 for verifing quantum-state 
steering. To illustrate what this condition implies in terms 
of the displacement sensitivity of such an optomechanical 
device, in the left panel of Fig. [2] we compare various dis- 
placement noises with respect to the Standard Quantum Limit 
(SQL)— sf QL = 2h/(mQ. 2 ) |26). We introduce three char- 
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acteristic frequencies of these noises where their displace- 
ment spectrums intersect the SQL as benchmark through — 
£lp = [2mK m k B T /{hm)] l l 2 , Cl q = [a 2 /(Tim)} 1 / 2 , and £2 A = 
Q q [2r\/{\ - tj)] 1 / 2 . In the right panel of Fig.|] we learn that 
the thermal noise from thermal fluctuation, and the sensing 
noise from optical loss and quantum inefficiency, need to be at 
least below the SQL in order to prepare and verify quantum- 
state steering. Note that the steerability does not depend on 
the absolute value of the noise spectrum; one can therefore 
have the flexibility to choose the appropriate frequency range 
to carry out the experiment, depending on the specific setup. 
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Notes. — While preparing this draft, an interesting protocol 
for realizing steering with pulsed optomechanical devices, in- 
stead of using a steady continuous optical field considered 
here, has been proposed by He and Reid l27l . 
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Appendix A: Wiener-Hopf method 

In this appendix, we show how to derive C T B C in 
Eq. (24) , which is equivalent to solving an integral equation 
shown in Eq. (28) , with the Wiener-Hopf method. We first 
show the general formalism and then specialize to the large- 
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bandwidth and strong-measurement limit that we have consid- 
ered. 



1. General formalism 

Here we first discuss the general formalism. In the 
continuous-time limit, C is a 2 x 2 matrix with elements given 
by Cij(t) and, similarly, the elements of B are Bij(t—t') = 
1,2) and they only depend on time difference due to station- 
arity). We define K = B _1 C or equivalently, BK = C, and 
K is a 2 x 2 matrix with the elements satisfying the following 
integral equations: 



2 A) 

£ / dt'B ik {t-t')K kj (t')=Qj(t), 

k=l J-°o 



(Al) 



and therefore, inverting B is equivalent to solving the above 
integral equation. Note that the upper limit of the integration 
is 0, instead of +°° in which case it can be solved simply by 
using Fourier transform. This arises naturally in the classical 
filtering problem with only past data that is available. The 
procedure for using Wiener-Hopf method to solve this set of 
integral equations goes as follows. 

Firstly, we extend the definition of Kij(t) and Cy(f) to t > 
but requiring Kij(t) = Cy(f) = if t > 0, namely 



Kij(t) -> Kij(t)&(t), Qj(t) -> Ctj(t)0(t) , 



(A2) 



This allows to extend the upper limit of the integral to be +°° 
without changing the result. 

Secondly, we apply the Fourier transform 



f{(0)=F\f{t)V 
of the above equation and obtain 



dte imt f{t), 



£B /(fe (©)^(0))-Qy(©) 

k=\ 



0. 



(A3) 



(A4) 



where [/(©)]_ means the part of f{(0) that is analytical (no 
poles) in the upper-half complex plane by using the following 
decomposition: 



f{(0)=\f{(0)} + + [/(©)]_, 



(A5) 



and [f(co)} + is the part that is analytical in the lower-half 
complex plane. From the definition of Fourier transform in 
Eq. ( |A3| l, the inverse Fourier transform of [/(o))]_ vanishes 
for / > from the residue theorem, namely 



[/(©)]_ =/(o@(0 = o, w>o. 



(A6) 



Let us focus on the equations associated with the first col- 
umn of C, i.e., 7 = 1 (the situation for j = 2 will be similar), 
and we rewrite Eq. ( |A4[ ) explicitly in terms of their compo- 
nents: 

[Bll(£»)£n(£») +Bi 2 (£»)£ 21 (a)) -Cn (fl))]_ =0, (A7) 
[B 2l (co)K n (co)+B 22 (co)K 2l (co)-C 21 (co)]_=0. (A8) 



Thirdly, if B\\ (t ) = B\\(—t), we can factorize Bn(ffl) as 

fin(fl>) = 0+(<D)0_(fl>) (A9) 

which is the Fourier counterpart of the Cholesky decomposi- 
tion in time domain. We now can express K\\ (co) in terms of 
K 2 i (to) in Eq. ( fA7) . We use the fact that 

[/»]_=() =► [/(£»)£+((»)]_ =0, V| (A10) 

Multiplying Eq. ( |A7| > by (p^ l (co), we get 



1 

9- 



Cn 

9+ 



Bj 2 K 2l 



(All) 



Plugging it into Eq. ( |A8) , we obtain 

B21 ' ' 



B 22 



B 2l B l2 
fill 



^21 + 



(P- 
_l_ BjA 



where we have used the fact that: 



Cn 
<P+ 

B i2 K 2l 



C 21 



= 
(A12) 

(A13) 



Again, if B 22 — B 2 \B\ 2 /B\\ is an even function of time — 
due to stationarity and time-reversal symmetry. We can make 
a similar factorization to the one shown in Eq. \A9\: 



(A14) 



B\ 2 K 2 \ 


B\ 2 K 2 \ 


+ 


B\ 2 K 2 \ 


<P+ 




+ 


9+ _ 



Buico) ~ y =y+(a>)y-(Q))- 

Bu(co) 



The same as the approach for deriving Eq. ( |A1 1[ ) by using the 
fact shown in Eq. ( |A10[ ), we obtain 



1 

•^21 — ^r- 

V- 



1 

¥+ 



C 2 i 



Bn 
<f>- 



Cn 
9+ 



B^ 
9- 



B\ 2 K 2 \ 



(A15) 

Note that in the above equation, K 2 i appears in both sides of 
the equation, which might seems to be difficult to solve. Ac- 
tually, since K 2 \ is analytical in the upper half complex plane, 
[Bi 2 K 2 \/(p + ] + only depends on the value of K 2 i at the poles 
of B12/9+ m me upper-half complex plane. One only need 
to solve a set of simple algebra equations by evaluating the 
above equation on these poles. 



2. Large-bandwidth and strong-measurement limit 

In this section, we consider the case of the large-bandwidth 
and strong-measurement limit, and, for the oscillator, we have 

x(t) = J'jt'GAt-Ol-aa^+F^t')}, (A16) 

p(t)=mf df'G. t (r-0[-aai(f')+^th(0]> (A17) 
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and, for the output optical field, we have 



where 



(A18) 



Here 



/I7 [«2 



(r) + (a/ft)jc(0] + v / l-i7«2(0- ( A19 > 



G,(f) = — e 



-Kmt/2 



sin <a,„f 



(A20) 



is the Green's function of the mechanical oscillator, and in the 
strong-measurement limit — the frequency at which we carry 
out the measurement is much higher than that of the mechan- 
ical frequency, the oscillator can be treated as a free mass and 

G.v (01 free mass =t/m. 

By using the fact that 

(dj(t)d k (t')) sym = {hj{t)h k {t')) sym = -8 jk S(t-t') (A21) 

for j, k = 1,2, and 

(F tb (t)F th (t')) sym = 2mK m k B T8(t - 1') , (A22) 

we obtains the elements for covariance matrix B of {y\ , yj) in 
the frequency domain (spectral density): 



Bn(0)) = l 
Bi 2 (ft>) = - 
B 2 i{(0) = - 



a 2 



a 2 ~ 



■q—G x (co)., 



a 2 



B 22 (co) = 1 + r)- T S xx ((o), 
h 

and the correlation between (y\ , y 2 ) and (i(0) , p(0)): 
Cu{(o) = -^r]aG* x ((o), 



(X 

C 2 i(ffl) = y/rj -S xx (co), 
n 

cc 

C 22 ((o) = imQ.y/rj-S xx (co) , 
n 



(A23) 
(A24) 

(A25) 

(A26) 



(A27) 
(A28) 

(A29) 
(A30) 



S xx {(0) = \G x {co)\ 2 (a 2 +4mK m k B T 



(A31) 



and the Fourier transform for the mechanical Green's func- 
tion is given by [strong-measurement limit is taken by setting 
7,„, C0 m -> 0]: 



G x {co) 



-1 



m(a> 2 - co 2 +iK m co) ' 



(A32) 



In this case, we can easily carry out the factorization. For 
q>±(oo) [cf. Eq. jA9\ ], we have 



m + (Q)) = ffl_(q)) = l, 
For y/±(co) [cf. Eq. flA14| i], we have 



(A33) 



a- 



y + ((o)y-((o) = l + Ti- T \G 2 (co)\ 2 Ul-ri)a 2 +4mK m k B T] 
^ + {Kl-2co 2 m )co 2 + coi + 2{a 2 /hm) 2 ^ F 



(0 4 + (K 2 -2(0 2 )(0 2 + (0i 



This leads to 



V+((o) = f-(co) = 



(o - a>i)(a> - ©2) 
((0-(a[)((0 ~(o' 2 )' 



(A34) 



(A35) 



where (Oj and co'j (j = 1,2) are the roots of the numerator and 
denominator of Eq. ( |A34| i in the upper-half complex plane, re- 
spectively. Given the expression for (p± and \jf±, we can solve 
Kij by using Eqs ( |A1 1 J > , and Eq. ( |A 1 5[ >. This in turn allows us 
to obtain V. s = A — C T K and in the time domain, it reads: 



(V s )y=Ay-£/ dt'C ki (t')K kj (t'), (A36) 



from which we obtain Eq. d32b. 



